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Abstract 

We present a first principles investigation of the lattice dynamics and 
electron-phonon coupling of the superconductor MgB2 and the isostructural 
AIB2 within the framework of density functional perturbation theory using 
a mixed-basis pseudopotential method. Complete phonon dispersion curves 
and Eliashberg functions a 2 F are calculated for both systems. We also report 
on Raman measurements, which support the theoretical findings. The cal- 
culated generalized density-of-states for MgE>2 is in excellent agreement with 
recent neutron-scattering experiments. The main differences in the calculated 
phonon spectra and a 2 F are related to high frequency in-plane boron vibra- 
tions. As compared to AIB2, they are strongly softened in MgB2 and exhibit 
an exceptionally strong coupling to electronic states at the Fermi energy. The 
total coupling constants are \MgB 2 = 0-73 and \aib 2 = 0.43. Implications for 
the superconducting transition temperature are briefly discussed. 
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Recently superconductivity has been observed in MgB2 with an exceptionally high tran- 
sition temperature for such a simple compound (T c ~39 K) [I]]. This system has in the 
mean time received a lot of attention from many experimentalists as well as theorists since 
understanding superconductivity in this binary compound should be much easier than in 
the high-T c cuprate materials studied extensively for more than ten years now. This is due 
to the simpler lattice structure and the missing complication due to magnetism and strong 
electron-electron correlations. Already early experiments a s well as theories f||| have 
suggested that we are dealing here with an s-wave superconductor based on strong electron- 
phonon coupling although alternative coupling mechanisms have been discussed too |jlQ| , [TTj| . 

Modern band-structure calculations together with ab initio determination of the phonon 
dispersion, electron-phonon coupling and solution of the Eliashberg equations to calculate 
the transition temperature should be the way to help understand the role played by the 
electron-phonon coupling in this material. Although during the past six weeks many con- 
tributions have been published concerning the electronic bandstructure P,|9|.|T2|j only a few 
attempts have been made to proceed along the above mentioned roadmap to a deeper under- 
standing of the electron-phonon coupling. Most investigations are restricted to calculations 
of the phonon modes for the T-point only P,|T5| . Estimates of the coupling strength are there- 
fore not very accurate thus asking for a more complete treatment of the electron-phonon 
coupling in the whole Brillouin zone (BZ). 

We have carried out such a systematic study of the lattice dynamics and the electron- 
phonon coupling using the mixed-basis pseudopotential method. To get a better understand- 
ing of the relevance of the electron-phonon mechanism we have studied two isostructural 
systems: MgB 2 , which is superconducting, and A1B 2 , for which no superconductivity has 
been found so far. This offers the possibility to compare the phonon dispersion curves to 
identify those which are strongly influenced by electron-phonon coupling even without cal- 
culating the coupling strength. Furthermore, the calculation of the Eliashberg function a 2 F 
for both systems allows for a consistency check of the proposed phonon-mediated pairing 
mechanism. 
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Only very recently after completion of our study we became aware of a similar investiga- 
tion by Y. Kong et al. [[TJ]] which calculated the full phonon dispersion, the electron-phonon 
coupling and T c for MgB2 within the LMTO-method. Our calculations are very similar to 



those presented in Ref. however, due to the parallel treatment of MgB2 and AIB2 they 
offer additional information. 

In contrast to all other calculations presented so far we have structurally optimized both 
systems, thus we remain fully in the framework set by the theory. For comparison and checks 
of sensitivity of our results we have also studied certain phonon modes using the experimen- 
tal structure. All our calculations are carried out using the mixed-basis pseudopotential 
method. For Mg and Al we have used well tested pseudopotentials of Martins- Troullier- 



15] and BHS type [16]. The boron pseudopotential was constructed according to the Van- 



derbild description [pj]] which led to a fairly deep p-potential for the boron atom which 
however could be dealt with very efficiently due to the mixed-basis formulation. The wave 
functions were constructed from localized s and p functions at the Boron sites supplemented 
by plane waves with an energy cut off of 16 Ry. Detailed tests were carried out to assure 
convergence with respect to the number of plane waves as well as with respect to k-point 



sampling flT8|]. Different Monkhorst-Pack k-point sets have been used up to 13824 k-points 



in the BZ together with a Gaussian smearing of 0.2 eV. The calculation of the phonon 
dispersion is based on a recently developed mixed-basis perturbation approach ]19|] which 
also allowed for the calculation of the electron-phonon coupling p0| |. Finally, estimates of 



T c are obtained by solving the linearized form of the Eliashberg equations ||21|| . Structural 



parameters were found to be fairly insensitive to the k-point sets while individual phonon 
modes were very sensitive to the sampling thus pointing to strong electron-phonon coupling. 
All our calculations were carried out using the local exchange-correlation potential of Hedin 
and Lundqvist |22]j . 

Since both systems crystallize in the so-called "AIB2 structure" with alternating hexago- 
nal (Mg,Al)-layers and graphite-like B-layers structural optimization required optimization 
of volume V and c/a ratio. Table 1 gives calculated results for V and c/a as well as the 



bulk modulus B together with the experimental values. For both systems a slightly smaller 
volume has been found by the theoretical treatment compared to the experimental one. 
This is a typical behavior in LDA-calculations. The c/a ratio is slightly larger. The elec- 
tronic bandstructure obtained for MgB2 is very similar to the ones obtained by other groups 
PPP^] with hole pockets around the T-point extending along the T-A-direction. A1B 2 has a 
very similar bandstructure, however, due to the additional electron density no hole pockets 
around T and the T-A-direction exist. 

In Table 2 we have summarized our phonon results for the T-point obtained from cal- 
culations based on the experimental structure as well as on the optimized one. With the 
exception of the E 2ff -mode (in-plane boron mode) all other modes are fairly insensitive to 
structural changes as can be seen from Table 2. The same holds also true for the sensitivity 
with respect to k-point sets. All previous calculations of the T-point phonon modes for 
MgB2 agree very well with these results with the exception of the E 29 -mode where shifts 
of ±100 cm -1 have been observed P Jl~^JT4]| . As extensive studies have shown this is 



mostly an effect of k-point sampling and indicates strong electron-phonon coupling for this 
mode. Comparing these results with those obtained for A1B 2 the most dramatic changes 
occur again for the E 2g mode which stiffens substantially and becomes the highest one, and 
the Big mode (out of plane boron mode) which is strongly softened. 

We have measured the E 2g mode by Raman scattering from samples of commercially 
available MgB 2 and A1B 2 powders. The samples showed clear x-ray diffraction patterns and 
consisted of crystalline grains up to 10 jttm (MgB 2 ) and 50 fim (A1B 2 ). Measured spectra 
showed only one prominent line in agreement with space group P6/mmm and were not par- 
ticularly polarization dependent. Typical results are presented in Fig. 1. The experimental 
result for A1B 2 is in almost perfect agreement with the calculated one. For MgB 2 the Raman 
line is very broad, however, the peak position agrees reasonably well with the calculated one 
while the linewidth is likely to be due to strong electron-phonon coupling (see discussion 
below). 

Having studied carefully the T-point modes we calculated also the full phonon dispersion 



for both systems. Results are presented in Fig. 2 which were obtained by determining the 
dynamical matrix on a (6,6,6) reciprocal lattice grid in the hexagonal Brillouin zone. The 



results for MgB2 agree very well with those of Ref. [I3[] and give thus additional credibility 
to the theoretical treatment. 

Comparing the phonon density-of-states shown in Fig. 2 certain features are worth men- 
tioning. The biggest difference between the two systems is showing up in the intermediate 
region of phonon frequencies where AIB2 shows nearly a gap due to the strong dispersion of 
the relevant modes while for MgB2 this gap has nearly completely disappeared. High density 
of states at the upper and lower end of the frequency range show up in both systems. To 



compare with recent measurements by inelastic neutron scattering [23] we have calculated 
also the generalized density-of-states (weighted by the inelastic scattering cross section and 
the mass) broadened by a Gaussian with a width of 4 meV. The result in Fig. 3 should be 



compared directly with Fig. 1 in Ref. [fj4j. An almost perfect agreement with respect to 
peak positions, shoulders and even with the relative ratio of contributions is obvious. We 
don't find any indication of a peak in the spectra near 17 meV, as observed by Sato et al. 
|25|. 



In the dispersion curves shown in Fig. 2, the layered structure of the crystals is reflected 
in a weak dispersion of the optical branches along T-A and in an anisotropy in the slopes 
of the acoustic branches in agreement with experiments |p6| . Besides many similarities two 



very significant differences between A1B 2 and MgB 2 can be seen. The first one is related to 
the branches which evolve from the doubly degenerate E 29 -mode at the T-point. In MgB 2 , 
these branches are strongly renormalized towards lower frequencies mostly near V and along 
the T-A direction. This is probably related to the hole pockets found in the electronic 
bandstructure of MgB 2 which are absent in A1B 2 . The second distinctive difference is the 
behavior of branches starting from the B lfl -mode which in MgB 2 are significantly harder in 
certain regions of the BZ than in A1B 2 . 

To address the superconducting properties we have calculated the so-called Eliashberg 
function a 2 F(u), using a very fine (36,36,36) k-point grid in performing the Fermi-surface 



average of the electron-phonon matrix elements p7| . Results are presented in Fig. 4. For 
MgB 2 we find indeed a very large contribution to a 2 F in the intermediate region between 
60 and 70 meV, which is mainly due to strongly softened in-plane vibrations of the boron 
atoms as mentioned above. At T, for the E2 S -mode this strong coupling results in a very 
large linewidth of 15 meV in accordance with the broad feature seen in the Raman spectrum 
(Fig. 1). For AIB2, in contrast, the biggest contributions show up at very high frequencies as 
well as in the regime of the acoustic modes, however, they are substantially smaller than the 
main contribution in MgB2. From the Eliashberg function we calculate the electron-phonon 
coupling constant 



a CO 



which gives \MgB 2 =0-73 and Aazb 2 =0.43. For the logarithmically averaged phonon frequen- 



cies as defined in Ref. [14] we find 60.9 meV and 49.9 meV for MgB2 and AIB2, respectively. 



The values for MgB 2 agree well with the results of Ref. []I3 ] . The fairly large value of A for 
A1B 2 however is unexpected. 



Within the dirty limit of superconductivity [F2SJ the calculation of the transition tem- 
perature T c requires the knowledge of a 2 F and the Coulomb pseudopotential fi* As 
commonly done, we treat \x* as an adjustable parameter. A T c of 40 K for MgB2 is obtained 
for /x*=0.05. The same /x*-value leads to T c ~ 9 K for the system A1B 2 . This result is in 
contrast to the experimental situation where no superconductivity has been found so far for 
A1B 2 . There are basically two possibilities to reconcile the theoretical results with the exper- 
imental situation. One possibility is a different screening in AIB2 compared to MgB2 which 
has been speculated about in Ref. p9| , leading to different values for \i* . Alternatively, the 



approximation of an isotropic superconductor (dirty limit) might not hold. Due to the fact 
that we have dealt with both systems on equal footing new questions have emerged which 
need further studies. 

We have presented here first principles calculations of the phonon dispersion and electron- 
phonon coupling for two systems MgB2 and AIB2 which crystallize in the same lattice struc- 
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ture, however, which have fairly different phonon dispersion curves. These results are in 
excellent agreement with measured quantities. The calculation of the superconducting tem- 
perature, however, has a problem which might be due to approximations involved or due to 
the fairly restricted knowledge about the Coulomb pseudopotential fi*. 

The authors would like to thank Dr. Reichardt, Dr. Pintschovius and Dr. Schweiss for 
many helpful discussions. 
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TABLES 

TABLE I. Structural parameters of the optimized geometries. Experimental values (in brack- 



ets) are taken from Refs. [ 30, 26| . 



v (A 3 



c/a 



B (Mbar) 



MgB 2 
A1B 2 



28.481 (29.010) 
24.617 (25.578) 



1.153 (1.142) 
1.09 (1.084) 



1.47 (1.20) 
1.84 



TABLE II. Comparison of the calculated L-point phonon frequencies for the experimental (I) 
and optimized (II) geometries. Values are given in meV (cm -1 ). 



Mode 


MgB 2 

I 


II 


I 


A1B 2 

II 


Eln 


39.9 (322) 


40.5 (327) 


33.0 (266) 


36.6 (295) 


A 2u 


48.9 (394) 


50.2 (405) 


48.6 (392) 


52.1 (420) 


E 2ff 


66.5 (536) 


70.8 (571) 


118.3 (954) 


125.0 (1008) 


B i<? 


86.3 (696) 


87.0 (702) 


59.9 (483) 


61.3 (494) 
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FIGURES 

FIG. 1. Micro-Raman spectra obtained from polycrystalline grains of MgE>2 and AIB2 at room 
temperature (A = 514.5 nm) 

FIG. 2. Theoretical phonon dispersion curves along high-symmetry lines of the hexagonal BZ 



(notation after |31|) and density-of-states (DOS) of MgB2 and AIB2. Dots represent actually 



calculated modes, lines are obtained by Fourier interpolation. 

FIG. 3. Calculated generalized phonon density-of-states of MgB2- Values for the incoherent 
scattering cross sections are taken from Ref. 



FIG. 4. Calculated Eliashberg functions a 2 F(u)) for MgB2 and AIB2. 
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